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The ideas of the hnear combination of atomic orbitals (LCAO) method, well known from the study 
of electrons, is extended to the classical wave case. The Mie resonances of the isolated scatterer in 
the classical wave case, are analogous to the localized eigenstates in the electronic case. The matrix 
elements of the two-dimensional tight-binding (TB) Hamiltonian are obtained by fitting to ab initio 
results. The transferability of the TB model is tested by reproducing accurately the band structure 
of different 2D lattices, with and without defects, thus proving that the obtained TB parameters 
can be used to study other properties of the photonic band gap materials. 



PACS numbers; 42.70Qs, 41.20Jb, 71.15-m 
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In recent years experimental and theoretical studies of 
artificially manufactured periodic dielectric media called 
photonic band gap (PBG) materials or photonic crys- 
tals, have attracted considerable attention PBG 
materials can have a profound impact in many areas in 
pure and applied physics. PBG materials are often con- 
sidered as analogous to electronic semiconductors. The 
existence or not of spectral gaps in periodic PBG ma- 
terials or localized states in disordered systems, in anal- 
ogy with what happens to the electronic materials, is 
of fundamental importance. While it is by now firmly 
established that certain periodic arrangement of dielec- 
tric structures has a full PBG in 2D and 3D it 
is not clear what mechanism is responsible for the for- 
mation of gaps. The relative importance of the roles of 
the two different mechanisms, single scatterer resonances 
and macroscopic Bragg-like resonances, in the formation 
of gaps and locahzed states is still being debated. Prelim- 
inary results have shown that there is a direct corre- 
spondence between the gaps calculated by plane wave ex- 
pansion and the Mie resonances Q of an isolated sphere. 
It is surprising that the positions of the Mie resonances 
approximately coincide with the center of the bands. It 
is tempting to suggest that the Mie resonances of an iso- 
lated scatterer to play the role of the energy levels of an 
isolated atom in a crystal. Also, it would be possible 
to formulate the problem in a simpler way, similar to the 
tight-binding (TB) formulation of the electronic problem. 

It is well known that the TB method has proven to be 
very useful in studying the electronic properties of solids 
[^-^. In an empirical TB approach, matrix elements 
of the Hamiltonian between orbitals centered on differ- 
ent sites are treated as parameters which are adjusted 
to obtain the band structure and the band gaps, which 
have been determined by other more accurate methods. 
The parameters obtained in this way are then used to 
study other properties of the systems, such as surface 
states, impurities and properties of disordered systems. 
The success of the TB formulation has been tested in the 
studies of Si, C and hydrogenated amorphous systems 



In this paper, we show that it is possible to extend the 
ideas of the LCAO method to the classical wave case. 
The Mie resonances of the isolated scatterer in the clas- 
sical case play the role of the localized eigenstates in the 
electronic case. However, there exist two important dif- 
ferences. First, Mie resonances' states are not localized, 
in fact they decay too slowly, as l/r as r — > oo and 
this may lead to divergences in some matrix elements. 
However, in a lattice environment they become localized, 
with a localization length analogous to the interparticle 
dimension. Second, in the classical wave case, as opposed 
to the electronic case, the host medium supports prop- 
agating solutions for every frequency. For large wave- 
lengths, this is the dominant propagation mode since no 
resonances have been excited yet, while for wavelengths 
comparable to the particle dimension, transmission is 
achieved through coupling with the fully excited local- 
ized resonances. In that sense, we may assume that a 
plane wave is coupled with the lowest frequency band, 
while all higher bands are described solely by the higher 
excited resonances. This assumption is mostly justified 
in the case of wide gaps and narrow bands, which is the 
one that we will study. We find that under appropriate 
rescalings, all matrix elements scale well with distance, 
and we obtain their functional dependence. This TB pa- 
rameterization of PBG materials will be very useful in 
studies of impurities and the effects of disorder in these 
materials. 

We will contend ourselves in the scalar wave equation 
case of a 2-D periodic array of N infinitely long dielec- 
tric cylinders, with periodic boundary conditions and the 
incident plane wave E-polarized. We assume the normal- 
ized electric field for each band to be given by 
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where n = 0, 1,2, ... is the resonance's (or band) index, 
^'„(r — R) stands for the wavefunction of the n'th res- 
onance localized at R, and has an angular symmetry 
^' ^ cos(7i6'). = 0, = 1 for 71 7^ and are functions 
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of the frequency (fc = |fc|) with |c^p + |c^p = 1, r, R, 
k are 2-D vectors, and we have assumed a unit area unit 
cell. We will assume a frequency independent normal- 
ized resonance function ^E", with J ^'*^(r)^'„(r — R)df = 

Smn^iR), as weh as / \E'*j(f— R)e^^^dr = so that 
J E^Endf = Smn in order to simplify the problem and 
make a better correspondence with the electronic case. 
This will turn to be a good approximation for our case 
of interest. For the lowest frequency band, we should ex- 
pect Cq — > 1, Cq — > for |A:| ^ and Cq — > 0, c§ 1 for 
\k\ - \G\/2. 

The "Hamiltonian" for the scalar wave equation is, H = 
— V^/e(r) and the eigenfrequencies liP' ji? of the system 
can be found by diagonalizing the Hamiltonian matrix 
ifmn = / E^^HEndr. Then : 

^oo-|4pi|^ + |cgp (^Poo + J2^oo{R)e^'^^ (2) 

where we have neglected crossing terms between the 
plane wave and the resonances (this is not completely 
justifiable, but it can be argued that they are not impor- 
tant for frequencies close to the gap, and so for impu- 
rity and disorder studies). Poo = f q {r ) H q {r ) dr a.nd 
700 (^) = J ^'*(r)7J4'o(f- R)df, and (e) is the averaged 
dielectric constant. We argue that the functional form 
of |co(/c)P is similar to the form of the scattering cross 
section of a single cylinder for the n = (or s-wave) case, 
so that |cJ(A:)P ~ e'^^^X. Here = \k\c/ {coo , luq 
is the single cylinder Mie resonance frequency, A(/) is a 
function of the fiUing ratio / of the form A(/) = bi/f''^. 
The power fj, has to be larger than 2 in order to preserve 
the correct slope at |fc| 0. For simplicity, we choose 
fi = 4. 

The second band (n = 1 or p-like) has a ^ ~ cos 9 
symmetry, and will consist of two linearly independent 
polarizations, Px and py. There will be crossing terms 
between these two polarizations, but not between s-like 
and p-like resonances. Using the standard notation 
for the matrix elements, and taking into account only first 
and second nearest neighbors, then the nonzero Hamilto- 
nian matrix elements for the case of a square lattice with 
lattice constant a is: 

Hss^\cl\^\k\y(e) + \cines + 2VBx 

X (cos -I- cos 0y ) -I- 4:Vj^}. COS (j)x COS (j)y\ (3) 

-ffpxPx = + "^Vpll COS (py + 21/W cos + 

+2 (y^ll + ) cos 0, cos 0, (4) 
Hp.Py = 2 (V^^l - V^^l) sin 0, sin (5) 

where (j)^ = k^a, 4>y = kya and |fc| — y'fc^ -|- fc^ ^ 

T^*-^-* refer to the 7 matrix elements for first and second 
nearest neighbors respectively, e refers to the /3, on-site or 



diagonal, matrix element. Hp p is similar to Hp^p^ with 
X < — > V and Hr, n = H* „ . In this work we consider 

^ yyyx PxPy 

only these two bands. 

We have fitted the V and e matrix elements for the 
E polarization case with cylinders of dielectric constant 
e = 100 in vacuum, to the band structure of five dif- 
ferent rectangular lattices with large/small axis ratios: 
1, 1.05, 1.1, 1.15, 1.2 as well as to a hexagonal lattice, 
for six different filling ratios: / = .1, .2, .3, .4, .5, .6. 
We used a large value for the dielectric constant to ensure 
that we obtain large gaps and narrow bands. The qual- 
ity of these fits can be seen in Fig. 1 where we plot the 
bands as found numerically by the plane wave expansion 
(PWE) method along with the TB fit, for a square and a 
hexagonal lattice for two filling ratios. The perfect qual- 
ity of the fits is an indication of the potential usefulness 
of the TB method. 

We next plot some of the fitted matrix elements. The 
square root of the diagonal Sp^ and Sp^ matrix elements 
are plotted (Fig. 2a) as a function of the filling ratio f, 
while the off diagonal T^pir matrix elements are plotted 
(Fig. 2c) as a function of the dimensionless separation 
distance dij = Tij /a^ where is the separation distance 
between cylinders i and j and a is the cylinders' radius. 
Notice how bad the scaling is, especially for the Ipp^ ele- 
ment. Apparently, a satisfactory TB description can not 
be achieved with this two-center approach, but rather, 
with the inclusion of the lattice environment j|] as well. 

The proposed simple rescaling function (-D°")i for the 
diagonal matrix element that describes the lattice envi- 
ronment of cylinder i, and takes into account the filling 
ratio and the different symmetries is of the form: 

1 ^ rcos^(n6>,j) 

where 9ij is the angle between the symmetry axis of the 
p resonance on cylinder i and the direction, n=0,l,.. 
for the s,p,.. resonances, and the sum runs over the near- 
est neighbors of cylinder i. The power on the angular 
function was chosen so that the Px and Py resonances 
in the hexagonal lattice be the same. The only choices 
were 2 and 4, and it was found that 2 gives better re- 
sults. Eq. 6 is similar to what was used in Ref. |^ for the 
atomic orbitals, except for two differences: (a) here we 
take into account the resonance's angular symmetry, and 
(b) the exponentially decaying part is missing, reflecting 
the non-localized character of the EM resonances. Fi- 
nally, T — [7r/(a^/)] is a renormalization that takes into 
account that different structures, with the same lattice 
constant a and cylinder radius a, have different filling 
ratios, r is normalized so that r = 1 for the rectangular 
structures and r = 3/4 for the hexagonal. We will use 
this parameter only for the diagonal matrix elements. 

For the periodic case, the function (-D^")i is the same 
for every j, and will characterize the corresponding reso- 
nance. We find the diagonal matrix element to scale as 
= af,Ha5'(£)°")-<-Fa^(D°")-< where = luqu/c 
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is the corresponding dimensionless Mie resonance fre- 
quency. In Fig. 2b we plot ^Je^ vs the environment 
function (1/_D°"). We can see that y/e^ now scales very 
well, having a larger value for increasing lattice density. 
The same dependence is found for all bands. 

In order to rescale the off diagonal V matrix element 
between two neighboring resonances i and j, we need 
contributions from the neighbors that are close to the 
line joining i and j. Contributions have to be projected 
on the fij direction for the s resonance, while for the p 
resonance we have to project on its symmetry axis. Only 
first nearest neighbors will contribute. At the end, we 
have to normalize with the sum of all projection weights. 

A simple formula that describes the environment of the 
n resonance on cylinder i, along the fy direction is: 

where I runs over j's nearest neighbors (including j). 0"i- 
is the angle between the fu and fij directions for the 
s resonance (n = 0), and for the p resonance {n — 1), 
it is the angle between the i'th resonance's symmetry 
axis and the fu direction . Both angles are taken to 
have a range from — 7r/2 to 7r/2 |jl^. Finally, if we 
include screening in our considerations, then the ac- 
tual matrix element V to be used in a particular prob- 
lem, can be obtained by the fully rescaled one V, by 
V'J = V'J(1 " S'^)/[{D°ff)-^ + {D°^^)~^] where S^i 
is the same screening function used in Ref. S^^ = 
tanh(^6iX;i5^»,je-''^[(''''+'*^')/'^'jl'''), and is different for 

different matrix elements. The fully rescaled matrix el- 
ements are found to scale with separation distance as 
V^j = cid'j''^ + C3^^,^■'=^ In Fig. 2d we plot the rescaled 
VppTT matrix element. We see now that it is a smouth 
function of separation distance, except for the second 
nearest matrix elements which do not scale very well for 
large filling ratios (small distances). Apparently a screen- 
ing function that depends on the filling ratio as well is 
needed, but for simplicity and to keep some accordance 
with the electronic models used so far, we will use this 
one. 

We find that all matrix elements now scale very well 
with separation distance, so a satisfactory TB description 
of PBG materials is in order. The rescaling parameters 
used and the parameters for the rescaled matrix element 
curves' generation are presented in Tables 1 and 2. We 
have also found = 1-65 for all n, and bi = .068, 62 = 
1.23 for A(/). 



TABLE I. The parameters for the e elements. 
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0.0804 


0.0460 


0.716 


-0.0121 


5.000 




0.2371 


0.0890 


1.640 


0.0020 


0.320 



TABLE II. The parameters for the V elements. 
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0.075 


0.0008 


10.0 


-0.108 


4.00 


-0.00096 


1.30 




0.100 


0.00008 


13.5 


0.425 


6.14 


0.0550 


3.22 




0.700 


0.0015 


10.0 


-0.076 


5.40 


-0.0044 


2.32 



Next, we test our results on a nonperiodic lattice. The 
original periodic lattice is square with filling ratio / = 0.2 
and lattice constant a. We will solve the 3x3 supercell 
problem. Each supercell consists of 9 cylinders, and the 
whole lattice is periodic in terms of this supercell. We 
keep all cylinders in the supercell in place, except the 
middle one, which we move from it's periodic position 
(0, 0), a/4 towards the left, then a/4 towards the bottom, 
and then return to the periodic position (0, 0) along the 
diagonal (see insert graph in Fig. 3). Solving the super- 
cell TB problem involves a 9x9 matrix diagonalization 
for the s band, and a 18x18 matrix diagonalization for 
the p band. We plot the 3 highest eigenvalues of the s 
band, and the 3 highest and 3 lowest eigenvalues of the p 
band. Luckily, for the 3 highest eigenfrequencies of the s 
band, we do not need to take into account the free plane 
wave, since for those frequencies the resonances are fully 
excited, blocking the background propagation mode. We 
plot the eigenfrequencies for ka = (0,1/3), and directly 
compare our results with the ones obtained numerically 
by the PWE method, for the same exactly system. It is 
worth noting that in the PWE method, 100 plane waves 
per rod are needed for an acceptable numerical conver- 
gence. This is a factor of 10'* more CPU time than the 
TB method. 

The agreement is excellent, so our TB parameteriza- 
tion works very well for the impurity case too. The pho- 
tonic TB model proves to be an excellent tool for under- 
standing and describing the band structure and proper- 
ties of light in PBG structures. We hope that the calcu- 
lated TB parameters will be used in studies of 2D PBG 
materials and can stimulate enough interest to extend 
these studies in 3D photonic crystals. 
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FIG. 1. The first two frequency bands for a square lattice 
with filling ratio / = .1 (a) and f = A (b), and for a hexagonal 
lattice with / = .1 (c) and f = A (d). Circles correspond to 
numerical results from the PWE method, while solid lines 
correspond to the TB fit. 

FIG. 2. Fitted TB parameters for the second (p-like) fre- 
quency band, (a) y/e^ vs f . (b) ^Je^ vs the rescaled envi- 
ronment function 1/Z)°". Circles and squares correspond to 
yj£p^ and respectively of a rectangular lattice with the 
big axis along x, and triangles to a hexagonal lattice, (c) 
VppTT VS the separation d. (d) The rescaled VppTr vs d. Circles, 
squares and diamonds correspond to a rectangular lattice's 
Vppn elements along the small axis, the large axis and the 
diagonal, and triangles to a hexagonal lattice. All matrix 
elements are expressed in the dimensionless units of (wa/c)^. 

FIG. 3. The three edge eigenfrequencies for each band for 
a 3 X 3 impurity system. All cylinders are identical, with the 
midlle one moving from the equilibrium position, as pointed 
ill the insert graph. The wave vector is always constant 
ka = (0,1/3). The first band gap extends approximately 
from uja/c ~ .53 to uia/c ~ .87, while the second starts at 
uia/c ~ 1.03 
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Lidorikis et al. Figure 1 




Lidorikis et al. Figure 2 
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